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We investigate tlie deconfining pliase transition in SU(3) pure gauge theory and in full QCD with two flavors of 
staggered fermions by means of a gauge invariant thermal partition functional. In the pure gauge case our finite 
size scaling analysis is in agreement with the well known weak first order phase transition. In the case of 2 flavors 
full QCD we find that the phase transition is consistent with weak first order, contrary to the expectation of a 
crossover for not too large quark masses. 



1. INTRODUCTION 

To detect the deconfinement phase transition 
in pure gauge theories the expectation value of 
the trace of the Polyakov loop is used as an or- 
der parameter. In presence of dynamical fermions 
the Polyakov loop ceases to be an order parame- 
ter since Z(N) symmetry is no longer a symmetry 
of the action. Alternatively deconfinement could 
be detected by looking at Abelian monopole con- 
densation which can be detected P by means of 
the vacuum expectation value of a magnetically 
charged operator whose expectation value is dif- 
ferent from zero in the confined phase and goes 
to zero at the deconfining phase transition. 

In a different way, abelian monopole condensa- 
tion can be signaled by the free energy needed to 
create an abelian monopole in the vacuum. This 
free energy can be computed |2I3| through a gauge 
invariant lattice effective action which in turn is 
defined at zero temperature by means of the lat- 
tice Schrodinger functional and at finite temper- 
ature by means of a thermal partition functional. 
Since abelian monopole condensation is not re- 
lated, as the Polyakov loop, to a symmetry of 
the lattice action, it could still be useful in or- 
der to detect a phase transition even if dynam- 
ical fermions are included into the action. The 
lattice thermal partition functional in presence of 



a static background field is defined |3I4| as 



VUe- 



Uk{Lt,S) = Uk{0.x) = Ul^'-(x) 



(1) 



T — 1/aLt is the physical temperature. T he 
spatial links belonging to the time slice xt — 
are constrained to the value of the external back- 
ground field, 



Uk{x%,=o^U'k^'ix), (fc = 1,2,3), 



(2) 



U^^^{x) being the lattice version of the external 
continuum gauge field, the temporal links are not 
constrained. At finite temperature the relevant 
quantity is the free energy functional: 
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Zt[0] 



(3) 



When including dynamical fermions, the thermal 
partition functional in presence of a static exter- 
nal background gauge field, Eq. iQJ , becomes: 



VUe 



detM. 



Uk(Lt/Z)=U^(0,x) = Ul^''(x) 



where Sw is the Wilson action and M is the 
fermionic matrix. Notice that the fcrmionic fields 
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are not constrained and the integration constraint 
is only relative to the gauge fields: this leads, as in 
the usual QCD partition function, to the appear- 
ance of the gauge invariant fermionic determinant 
after integration on the fermionic fields. 

To detect monopole condensation we con- 
sider |2l5j the following quantity defined in terms 
of free energy F needed to create a monopole in 
the vacuum: 
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Zt[0] 
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In presence of monopole condensation F is finite 
and e~^/^p''y= 7^ 0. In practice it is easier to 
compute F'{(3) the /3-derivative of the free en- 
ergy (eventually, since F — at (3 = 0, F can be 
obtained by numerical integration of F'{(3)). 

2. SU(3) PURE GAUGE 

We consider SU(3) pure gauge in an abehan 
monopole background field. The gauge potential 
in the continuum is 



gb^x) 
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where n is the direction of the Dirac string and 
'^mon (integer) is the number of monopolcs. On 
the lattice the spatial links exiting from the 
sites at the boundary of the time slice xt = 
are constrained to the (lattice version of) 
the gauge potential Eq. jSJ. Links at spatial 
boundaries {xt ^ 0) can be constrained to the 
monopole background field ( "spatial fixed bound- 
ary conditions" ) or can be periodic ( "spatial peri- 
odic boundary conditions"). These two different 
choices (for the monopole field which vanishes at 
infinity) are equivalent in the thermodynamical 
limit. The simulations were performed on lat- 
tices of different spatial sizes (16'^, 24^^, and 32'^) 
and fixed the temporal extent {Lt ~ 4), using 
APEmille/crate in Bari. We find that our data 
for F'{l3,Lf^) can be fitted according to the fol- 
lowing scaling law 



F'{f3,Lf)^ 



(Lf)i/-(/3-/3,)-di 



(6) 



where the above scaling relation holds quite well 
for a large range of the scaling variable x = 
{Lf)^/''ij3 - Pc), with u = 0.334 ± 0.021 con- 
sistent with a first order phase (see Fig. 1). 
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Figure 1. We plot together data for F' (fbcs: 
"spatial fixed boundary conditions" , = Ls — 
2) and (pbcs: "spatial periodic boundary condi- 
tions", Lf = Ls). 

3. QCD {Nf = 2) 

We study QCD with two dynamical staggered 
fermions. We simulate the theory with a "cold" 
time-slice where (as in the case of pure gauge) 
the spatial links are constrained to the abelian 
monopole background field. In Fig. Q the 
derivative of the free energy with respect to the 
gauge couphng F'{P) is displayed together with 
chiral condensate data suggest that the peak in 
F'{(3) corresponds to the drop of the chiral con- 
densate. Using our data for F'{f3) on different 
spatial volumes and different bare quark masses 
we try to infer the critical behavior of two flavors 
full QCD near the phase transition. We varied the 
lattice size Lg {Lg = 16, 20, 32) and the staggered 
quark mass {rrig = 0.075,0.2676,0.5003). At 
fixed niq, as in the quenched case, the peak in- 




Figure 2. F' on a 16^ x 4 lattice in the case of 
full QCD with 2 dynamical flavors, is displayed 
together with chiral condensate. 

creases as the lattice volume is increased. At fixed 
spatial volume the critical coupling depends on 
ruq. All our lattice data can be described by a fi- 
nite size scaling function where the gauge critical 
coupling now does depend on the quark mass and 
is determined by the chiral critical point [SI7I5) 

F\P,Lf,m,) = 

ai{LfV (7) 

\{LfY/-{l3 ~ l3c{Q) - crr4) ~ di\" ' 

the relevant chiral critical exponents being com- 
patible with those of the three-dimensional 0(4) 
symmetric spin models |7] . We find that the crit- 
ical exponent v = 0.31±0.03 is consistent with a 
first order phase transition, even though weaker 
than in SU(3) pure gauge (analogous indications 
have been obtained by the Pisa group [HI)- 

4. CONCLUSIONS 

We investigated the phase transition in SU(3) 
pure gauge theory and in full QCD with two fla- 
vors of degenerate staggered fermions. To locate 
the phase transition we looked at the free energy 
to create a monopole in the vacuum. In the pure 
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Figure 3. F' {(3 , Lf^ , mq) rescaled by the factor 
{Lf^)'* . The values of Lf^ and iriq = am are 
displayed in the legend. 

gauge case our finite size scaling analysis is con- 
sistent with the well known weak first order phase 
transition. In the case of QCD we find that the 
phase transition is compatible with a weak first 
order phase transition (contrary to the expecta- 
tion of a crossover for not too large quark masses). 
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